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Field theoretical scheme of regular Big Bang in 4-dimensional physical space-time, built in the 
framework of gauge approach to gravitation, is discussed. Regular bouncing character of homoge- 
neous isotropic cosmological models is ensured by gravitational repulsion effect at extreme conditions 
without quantum gravitational corrections. The most general properties of regular inflationary cos- 
mological models are examined. Developing theory is valid, if energy density of gravitating matter 
is positive and energy dominance condition is fulfilled. 
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I. INTRODUCTION 

The Hot Bing Bang scenario built on a base of 
Friedmannian homogeneous isotropic cosmological mod- 
els (HICM) of general relativity theory (GR) is the foun- 
dation of relativistic cosmology. Well known physical 
laws allow to describe physical processes in evolving Uni- 
verse beginning from the time t\ ~ 10~ 4 s from the Bing 
Bang corresponding to the start of cosmological expan- 
sion (in the case of Friedmannian HICM the moment 
t = corresponds to cosmological singularity), when 
the energy density of gravitating matter was comparable 
with nuclear density and the temperature T ~ lO 11 ^. 
Observable confirmation of theoretical predictions of pri- 
mordial nucleosynthesis of light elements, explanation of 
large scale structure formation of the Universe are impor- 
tant achievements of the Hot Bing Bang theory. To de- 
scribe gravitating matter and physical processes in evolv- 
ing Universe at t < t\ one has to use modern theory 
of fundamental physical interactions, unified models of 
particle physics, that is the object of cosmology of early 
Universe. The creation of inflationary paradigm, which 
permits to solve a number of problems of standard cos- 
mological scenario, in particular, to explain the homo- 
geneity and isotropy of the Universe at initial stages of 
cosmological expansion is important achievement of early 
Universe cosmology 0. The further progress in cosmol- 
ogy of early Universe is connected with further perfection 
of particle physics theory and also theory of gravitational 
interaction. 

As it is known, GR has problems by description of 
gravitating matter at extreme conditions with extremely 
high energy densities, pressures, temperatures (problem 
of gravitational singularities). The Friedmannian HICM 
possess singular state in the past, in which some impor- 
tant physical characteristics of gravitating matter (en- 
ergy density, temperature etc) diverge and the scale fac- 
tor R(t) of Robertson- Walker metrics vanishes. As a re- 
sult the existence of Friedmannian HICM is limited in the 
time in the past. Corresponding problem - the problem of 
cosmological singularity (PCS) - is one of the most prin- 
cipal cosmological problems, which does not have general 
solution still. From physical point of view, the PCS is 



connected with the fact that gravitational interaction in 
the case of usual gravitating systems with positive val- 
ues of energy density p and pressure p in the frame of 
GR as well as Newtons gravitation theory has the char- 
acter of attraction, but not repulsion. Although in the 
case of gravitating systems with negative pressure (for 
example, massive and nonlinear scalar fields) the gravi- 
tational interaction in GR can have the repulsion charac- 
ter, however, the PCS can not be solved bytaking such 
systems into account in the frame of GR . There 
were many attempts to resolve PCS in the frame of GR 
as well as other classical theories of gravitation (see 0, Q 
and Refs given herein) . A number of regular cosmological 
solutions was obtained in the frame of metric theories of 
gravitation and also other theories, in the frame of which 
gravitation is described by using more general geometry 
than the Riemannian one. In connection with this, note 
that the resolution of PCS means not only obtaining reg- 
ular cosmological solutions, but also excluding singular 
solutions of cosmological equations, as a result generic 
feature of cosmological solutions has to be its regular 
character. Moreover, gravitation theory and cosmologi- 
cal equations have to satisfy the correspondence principle 
with GR in the case of usual gravitating systems with 
sufficiently small energy densities and weak gravitational 
fields excluding nonphysical solutions. The greatest part 
of existent attempts to resolve the PCS does not satisfy 
indicated conditions 

Usually the Universe evolution in the frame of standard 
cosmological scenario is considered beginning from the 
time t greater than the Planckian time t p ~ lCP 43 s,when 
the energy density p and temperature T were smaller 



than the Planckian ones: 
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10 A'.Then according to general opinion quantum grav- 
itational effects are negligible and gravitational field can 
be described classically. On the contrary, at Planckian 
conditions (p > p p , T > T p ) one means that quantum 
gravitational effects must be essential, and possibly the 
PCS can be solved by taking these effects into account. 
So, in order to solve the PCS the idea of quantum birth 
of the Universe was introduced 0, @- There are sev- 
eral realizations of this idea in the frame of quantum 
cosmology by using Wheeler-DeWitt equation (see re- 
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view [8j). The appearing closed micro- universe is trans- 
formed then into macro-universe by virtue of inflation. 
Although singular state with divergent energy density is 
absent in such model, the problem of the beginning in 
the time of the Universe remains. The theory of quan- 
tum micro-universe was developed further in the frame 
of loop quantum gravity. The space quantization permits 
to continue solution for micro-universe of closed and flat 
type to the past and obtain bouncing solution containing 
the compression stage before the expansion @. Note, 
this scenario has some vaguenesses, in particular: if at 
compression stage the Universe was a macro-object filled 
by gravitating matter what physical factors lead to trans- 
formation of macro-universe into micro-universe before a 
bounce, what is a flat model with finite volume ? 

The PCS was discussed also in the frame of other 
candidate for quantum gravity theory - string theory/ 
M-theory El [H El El El E3 Some cosmological 
scenarios with evolution stage before the Big Bang - 
Pre-Big Bang scenario E3 an d ekpyrotic scenario E3 
- were proposed. Note that HICM in a four-dimensional 
Einsteinian frame, generally speaking, are singular |ld| . 
The obtaining of regular bouncing solutions is connected 
with violation of condition of energy density positivity 
for gravitating matter . So, in Ref. E3 non-local neg- 
ative scalar field potential is introduced, in Ref. [l5j 
a component with negative energy density is used. To 
build oscillating cosmological model, the specific neg ative 
scalar field potential was introduced in Refs. 0,0]- Al- 
though the energy density does not diverge by transition 
from compression to expansion, the singularity connected 
with vanishing of the scale factor R(t) remains. 

As it was shown in a number of our papers (see ja| and 
Refs given herein), the applying of gauge approach to 
gravitation permits to build regular Big Bang theory by 
satisfying positivity property of energy density for gravi- 
tating matter. The present paper is devoted to consider- 
ation of principal features of developing approach and to 
further analysis of regular inflationary cosmology. In Sec- 
tion 2 generalized cosmological Friedmann equations for 
HICM in the frame of gauge approach to gravitation are 
introduced. In Section 3 these equations are applied to 
analyze inflationary cosmological models filled by scalar 
fields and usual gravitating matter. As illustration of reg- 
ular Big Bang theory in Section 4 particular inflationary 
cosmological model is discussed. 



II. GAUGE APPROACH TO GRAVITATION AND 
GENERALIZED COSMOLOGICAL FRIEDMANN 
EQUATIONS 

As it is known, the local gauge invariance principle is 
the basis of modern theory of fundamental physical in- 
teractions. The theory of electro-week interaction, quan- 
tum chromodynamics, Grand Unified models of particle 



physics were built by using this principle. From phys- 
ical point of view, the local gauge invariance principle 
establishes the correspondence between certain impor- 
tant conserving physical quantities, connected according 
to the Noether's theorem with some symmetries groups, 
and fundamental physical fields, which have as a source 
corresponding physical quantities and play the role of car- 
riers of fundamental physical interactions. The applying 
of this principle to gravitational interaction leads, gen- 
erally speaking, to generalization of Einsteinian theory 
of gravitation. Metric theories of gravitation including 
GR, in the frame of which the energy-momentum ten- 
sor is a source of gravitational field, can be introduced 
in the frame of gauge approach by using the localization 
of 4-parametric translation group E3 [23 • Because the 
localized translation group is, in fact, the group of gen- 
eral coordinate transformations, from this point of view 
the general covariance of GR plays dynamical role. At 
the same time the Lorentz group (group of tetrad Lorentz 
transformations) in GR and other metric theories of grav- 
itation does not play any dynamical role from the point of 
view of gauge approach, because corresponding Noethers 
invariant in these theories is identically equal to zero |17j . 
If one means that the Lorentz group plays the dynam- 
ical role in gauge field theory and the Lorentz gauge 
field exists in the nature, in this case we obtain with 
necessity the gravitation theor y in the Riemann-Cartan 
space-time (see, for example, |l8l Ejj ). Corresponding 
theory is known as Poincare gauge theory of gravitation 
(PGTG). Gravitational field variables in PGTG are the 
tetrad (translational gauge field) and Lorentz connection 
(Lorentz gauge field); corresponding field strengths are 
torsion and curvature tensors. As sources of gravitational 
field in PGTG are covariant generalizations of energy- 
momentum and spin tensors. Unlike gauge Yang-Mills 
fields, for which the Lagrangian is quadratic in the gauge 
field strengths, gravitational Lagrangian of PGTG can 
include also linear in curvature term (scalar curvature), 
which is necessary to satisfy the correspondence principle 
with GR. 

The first attempt to apply the simplest PGTG - 
Einstein-Cartan theory - in order to solve the PCS was 
made in Refs. jHHj. By using some classical model 
for spinning matter, non-singular cosmological solutions 
were obtained. However, it was shown later, these solu- 
tions have model character and critically depend on spin- 
ning matter description; by another spin description (for 
example, by means of Dirac field) the cosmological singu- 
larity does not disappear. Moreover, because in the case 
of spinless matter the Einstein-Cartan theory is identi- 
cal to GR, all singular solutions for Friedmannian HICM 
are exact solutions of Einstein-Cartan theory of gravita- 
tion. The next step to apply the PGTG in order to solve 
the PCS was made in Ref.]^. In the frame of PGTG 
based on general expression of gravitational Lagrangian 
Lq including both a scalar curvature and different in- 



3 



variants quadratic in the curvature and torsion tensors 
gravitational equations for HICM were deduced, these 
equations lead to the followin g: ge neralized cosmological 
Fricdmann equations (GCFE) [34j: 



k_ 
R 2 



— In 

dt 



Ry/\l + a(p-3p)\ 



8n p+%(p-3pY 



3M p 2 l + a(p-3p) ' 



(1) 



R + R[lny/\l + a(p-3p)\ 



R 

4tt p + 3p - § (p - 3pf 



3M p 2 1 + a (p - 3p) 



(2) 



where k = +1,0,-1 for closed, flat, open models re- 
spectively, a is indefinite parameter with inverse dimen- 
sion of energy density, M p is Planckian mass, a dot de- 
notes differentiation with respect to time. (The system of 
units with H = c = 1 is used). Note that Eqs. 0-0 are 
valid also in the frame of the most gener al g auge theory 
of gravitation - metric-affine theory [3, l24j . From Eqs. 
0-0 follows the conservation law in usual form 



p + 3H(p+p)=0, 



(3) 



where H = is the Hubble parameter. Besides cos- 
mological equations 0-© gravitational equations of 
PGTG lead to the following relation for the torsion func- 
tion S 



4 dt 1 



a(p-3p)\ 



(4) 



Before we start discussing inflationary cosmology on a 
base of GCFE 0-0, let us make several remarks about 
these equations. At a — > the GCFE 0-0 are trans- 
formed into Friedmann cosmological equations of GR. 
The value of a -1 determines the scale of extremely high 
energy densities. Solutions of GCFE 0-© coincide 
practically with corresponding solutions of GR if the en- 
ergy density is small \a(p — 3p)\ <§; 1 (p ^ |p). The 
difference between GR and PGTG can be essential at 
extremely high energy densities \a(fi— 3p)\ > 1. Ul- 
trarelativistic matter (p = ^p) and gravitating vacuum 
(p = —p) with constant energy density are two excep- 
tional systems because Eqs. 0-© are identical to 
Friedmann cosmological equations of GR in these cases 
independently on values of energy density and 5 = 0. 
The behaviour of solutions of Eqs. depends es- 
sentially on the following restriction on equation of state 
of gravitating matter at extreme conditions: p > |p or 
p < ^p. Because of strong nucleon interaction we have 
for nuclear matter p > ip |25j : we will put that this re- 
striction is valid also for gravitating matter at extreme 



conditions, and the scale of extremely high energy densi- 
ties determined by a -1 surpasses nuclear density. From 
physical point of view, we can suppose that the value of 
a -1 is smaller than the Planckian energy density, but it 
will be shown later, the behaviour of solutions of GCFE 
does not depend on this assumption. 

The GCFE lead to restrictions on admissible values of 
energy density. In fact, if energy density p is positive 
and a > 0, from Eq.0 in the case k — +1, follows the 
relation: 



Z = 1 + a (p - 3p) > 0. 



(5) 



The condition is valid not only for closed and flat 
models, but also for cosmological models of open type 
(k = —1) (see below). In the case of models filled by 
gravitating matter with equation of state p = p(p), it 
is easy to obtain the solution of the system of Eqs. 
and in quadratures |26| ; these solutions are regular in 
metrics, Hubble parameter, its time derivative and have 
bouncing character. The transition from compression to 
expansion takes place by reaching limiting energy den- 
sity defined by the following condition Z — 0. At first 
the conclusion on possible existence of limiting energy 
density for gravitating systems, close by which the grav- 
itational interaction has the character of repulsion, was 
obtained in Ref.j2^]. Later the idea on limiting energy 
density as "universal law of the nature" was discussed 
by M.A. Markov [2^, and the value of limiting energy 
density was postulated to be equal to the Planckian one. 
According to the value of limiting energy density in 
our case depends on parameter a and equation of state 
for gravitating matter at extreme conditions and can be 
essentially smaller than the Planckian one . 



III. INFLATIONARY COSMOLOGICAL MODELS 
AND ITS PROPERTIES 

Now by using GCFE 0-0 we will study homoge- 
neous isotropic models filled by interacting scalar field <j) 
minimally coupled with gravitation and gravitating mat- 
ter with equation of state in general form p rn — p m {Pm)- 
(The generalization for the case with several scalar fields 
can be made directly). Then the energy density p and 
pressure p take the form 



p = -4> 2 + v + Pm (p>o), 



P=\4> 2 -V + p m , (6) 



where scalar field potential V — V{<j>, p m ) includes the in- 
teraction between scalar field and gravitating matter. In 
the most important particular case of radiation (ultra- 
relativistic matter) the expressions of V{4>,p m ) can be 
obtained by taking into account temperature corrections 
for given scalar field potentials and the following re- 
lation for energy density p m ~ T A . Because the form of 
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equation of state for gravitating matter changes by cos- 
mological evolution and we do not know explicit form of 
scalar field potential, our analysis will be made without 
its concretization. Our main aim will be to investigate 
properties of inflationary cosmological solutions for early 
Universe. By using scalar field equation in homogeneous 
isotropic space 

dV 

^ + 3H0 = -— (7) 



we obtain from Eqs. ®, ©, Q the conservation law 
for gravitating matter 

( dV\ 

Pm 1 + "5 + 3if (p m +Pm) = 0. (8) 

\ & Pm J 

By using Eqs.®-© the GCFE @J-© can be trans- 
formed to the following form 



HIZ + 3acj) 2 



3a 



Y 



3a- 



dV 



k o 8n 
Z 2 = 



R 2 



3M 2 



1 



1 



Pm + ^ + V + -a ( AV - <j) 2 + p m - 3p 



Z, (9) 



H 



Z + 3a 



Z + H Z 



l 5a 02 + ^ y 



9a Y 



1 + 



dp 



rn , Pm 



2 rfp m 1 + |£ dp 2 



p m <9 2 1M 9a (p m +p m ) 2 d 2 
' 2 375" ^ m " 4F ) 



\ dp 2 



!*«- («£ 2 + Iy) I -12ai/0|3a|^ ^ 2 + V 



5U 9 p m +p 
90 8 



1 



d 2 V 
d(pdp„ 



3a 



d 2 V „ fdV 



8tt 
3M| 



3 d/Om 

Z- 18a 
1 



dV 



If ■ \ 2 

(p m + 3p m ) + -a (^4y - 2 + p m - 3p, 



Z, (10) 



where Z 



1 + a (AV - </> 2 + p m - 3p m ) 



and y = 



^ ( 3 ^™ ~ 4y ) ~ x ) ■ The formula takes 



the form 



S 



3a 
2Z 



dV 



(11) 



The relation JSJ) determines the following restriction on 
admissible values of variables for scalar field and gravi- 
tating matter: 



\ 2 < AV + a" 1 



3p„ 



(12) 



Now let us introduce the 3-dimensional space P with axes 
(</>, <fi, p m )- The domain of admissible values of scalar field 



</>, time derivative <fi and energy density p m in space P 
determined by ifT^l) is limited by bound L defined as 

Z = or <j> = ± (AV + a- 1 + Pm ~ 3 Pm y . (13) 

From Eq. @ the Hubble parameter on the bound L is 
equal to 



dV ■ 



(14) 



-Y 



Let us consider the most important general proper- 
ties of cosmological solutions of GCFE lj§|l- ljTU|) . At first 
note, by given initial conditions for variables (</>, 4>,p m ) 
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and also in the case k = ±1 for R there are two differ- 
ent solutions corresponding to two values of the Hubble 
parameter following from Eq. 10: 



dV ■ 
-3a — 6 ± VD 



Z + 3a[ 



o*1 



(15) 



where 



D = 



3M2 



[Pn 



2*' 



V 



r ( 



W - (f> 2 + p m - 3 P „ 



2 ]Z-^ 2 Z\ (16) 



Obviously, the expression l|15fl for H± will be regular, 
if Y > 0. This relation is valid, in particular, for all 
models filled by gravitating matter with p m > and 
scalar fields with potentials applying in chaotic inflation 
theory. Unlike GR, the values of H + and H_ in GTG are 
sign-variable and, hence, both solutions corresponding to 
H+ and H- can describe the expansion as well as the 
compression in dependence on its sign. Below we will 
call solutions of GCFE corresponding to H + and as 
i/+-solutions and //--solutions respectively. Note that 
at asymptotics like GR the sign of 7J_ is negative and the 
sign of H + is positive. In points of bound L we have D = 
0, H + = H_ and the Hubble parameter is determined 
by ifTlj l. If initial conditions correspond to asymptotics 
of i/_-solution, then unlike GR by compression stage the 
derivative <fi does not diverge and by reaching the bound 
L the transition from //--solution to iJ+-solution takes 
place. In fact, by using the following formula for H±- 
solutions 



dV 



>Z±VD 



Z = 6a ■ 



Z + 3a 



-Y 



(17) 



it is easy to show that 



lim H ± =H L + 

z^o 



2ir p n 



3aM 2 



Ay 

2 1 



(18) 



From (|18|l follows that in points of bound L the deriva- 
tives H + and H- are equal and its values do not depend 
on the model type, as a result we have the smooth tran- 
sition from i/_-solution to //.^-solution on bound L, and 
corresponding cosmological solutions for all types mod- 
els are regular in metrics, Hubble parameter and its time 
derivative. Note, that in points of bound L conditions 
of uniqueness of solutions of Eqs. l(5)l-l|TU)l are not ful- 
filled, as a result there are specific solutions 0, , tra- 
jectories of which are situated on the bound L and have 
with H± -solutions common points, -solutions reach 
the bound L and i/ + -solutions originate from them, the 
surface Z = is envelope in space P for cosmological so- 
lutions According to Eq. ([TT| and |gj the function S 
has the following asymptotics for H + - and //--solutions 
at Z -> 0: 



1 Z 
lim S = — lim — ~ =pZ 

z^o 4 z^o Z 



(19) 



Unlike flat and open models, for which H + = H_ only in 
points of bound L and regular inflationary models include 
H + - and f/_ -solutions reaching bound L, in the case of 
closed models the regular transition from H_ -solution to 
i/+-solution is possible without reaching the bound L. It 
is because by certain value of R according to l|15fl H6|) 
we have H + = H_ in the case Z ^ 0. Such models are 
regular also in torsion . Regular inflationary solution of 
such type was considered in Ref . |29j . 



All discussed cosmological solutions have bouncing 
character, but in the presence of scalar fields a bounce 
takes place not in points of bound L [Z = 0). In or- 
der to study the behaviour of cosmological models at a 
bounce, let us analyze extreme points for the scale fac- 
tor R(t): Rq = R(0), H = H(0) = 0. (This means that 
in the case of H + -solutions H + o = and in the case of 
//--solutions -f/_o = 0). Denoting values of quantities 
at t = by means of index "0" , we obtain from (|9"|)- (|10fl : 



152^o 2 + 9a 2 — U 



^ ) 



8tt 
3A/f 



l- 2 1 / - 2 

Pmo + -j9o + Vo + -^a [AVq - 4> + p m0 - 3p m o 



Zq, 



(20) 



8tt 
3M% 



V 



1 1 / \ 2 

~ 2 ( Pm0 + 3p " l °) + 1 a ( 4V ° ~ ^0 + PmO - 3p m0 



-3a 



d 2 V 



dV\ 2 



Zq + 3a 



-Yo 



, (21) 
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where Z$ = 1 + a (4Vo — 0q + p m o — 3p m o) • A bounce point is described by Eq. (|2(jp. if the value of i?o is positive. 
By using Eq. (12011 we can rewrite the expression of Hq in the form 



-3a 



Ml 
d 2 V 



1 1 / \ 2 

o) + j«(- - 2 



Vo + g (PmO - P 



dV\ 2 



2k 



- Z } {Z + 3a 



(22) 



We see from (1221 unlike GR the presence of gravitating 
matter satisfying the energy dominance condition (p™< 
p m ) does not prevent from the bounce realization|37j. 
Eq. (|2t)|l determines in space P extremum surfaces de- 
pending on the value of a and in the case of closed and 
open models also parametrically on the scale factor Rq. 
In the case of various scalar field potentials applying in 
inflationary cosmology the value of -H+o or ii_o is posi- 
tive on the greatest part of extremum surfaces, which can 
be called " bounce surfaces" [3^ . By giving concrete form 
of potential V and choosing values of Rq, (f>o, 4>o and p m o 
at a bounce, we can obtain numerically particular bounc- 
ing solutions of GCFE for various values of parameter a. 
Like GR, if initial value of scalar field at the beginning of 
cosmological expansion is not small (</> > 1M P ) cor- 
responding solution is inflationary cosmological solution 
containing in addition to inflationary stage also compres- 
sion stage, transition stage from compression to expan- 
sion and a stage after inflation with oscillating regime for 
scalar field. For given scalar field potential properties of 
regular inflationary cosmological solutions depend on ini- 
tial conditions at a bounce and parameter a. Numerical 
analysis of such solutions in the case of simplest scalar 



field potentials applying in chaotic inflation was carried 
out in Ref. j^. Note that by equal initial conditions 
characteristics of inflationary stage in developing thery 
coincide with that of GR |32j . 

The analysis of GCFE shows, that some properties of 
cosmological solutions depend essentially on parameter 
a, i.e. on the scale of extremely high energy densities. 
From physical point of view interesting results can be ob- 
tained, if the value of a" 1 is much less than the Planckian 
energy density [i^, i.e. in the case of large values of pa- 
rameter a (by imposing M p = 1). In order to investigate 
cosmological solutions at the beginning of cosmological 
expansion in this case, let us consider the GCFE by sup- 
posing that 



a (w - 4> 2 + p m - 3p m ^j 



»1, 



p m + ^<p 2 + V^a(iV-<j> 2 +p m - 3p m ) 2 . (23) 



Note that the second condition (12: il) does not exclude that 
ultrarelativistic matter energy density can dominate at a 
bounce. We obtain: 



R 2 



4V 



H 



W + 2cj) 2 



3p„ 




2vr 
3M 2 



\AV - 



Pm ~ 3p 

77 



(24) 



H [W + 2cj) 2 + Pm - 3 P , 



Y ] \4V — 



Pm - iPm + H 



W - Wcf> 2 + Pm - 3 Prl 



3 y 9 Y 
2 21 + C 



2 J2 



dp m . Pm+Prn 8V\ 9 (p m + Pm) d 



dp n 



i 



dv \ 2 dp 2 



(3 Pm - 4V) 



x (W -4> 2 + Pm - 3 Pm ) - 18 (V + h^j 
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12Hcf)- 



dV , 9 Pm +p rn d 2 V Aid 

1 + IP™ + 21/J 



2 <90<9p„ 



3 d/9„ 



x (W - 2 + p m - 3p m ) + 3^ f 2 + \y 



d 2 v - 2 /av 

a^ - 



4V - + p m - 3p m 



+ Pm - 3p m . (25) 



According to Eq. (|24|l the Hubble parameter in considered approximation is equal to 



H+ = 



W + 2(f, 2 + Pm - 3 Pm + -Y 



4V- 



ftn - 3p 



r 



\/ 3M P 2 ( 



4^ - 2 + p m - 3p m ) - — 



(26) 



and extreme points of the scale factor are determined 
by the following condition 



i+9 
H 



2tt 



3M 2 



4Vq - O + Pmo - 3p m0 



4Vq - + PmO - 3p m0 



(27) 



From Eq. IT25H the time derivative of the Hubble param- 
eter at extreme points is 



f 2tt 



= | \ iVo - O + PmO ~ ZPmO 



-3 



d 2 V 



,2 / *0 



dV\ 2 



18 



(S)> 2 



4Vq - 00 + PmO - 3pm0 



W Q + 20 o + p m0 - 3pm0 + -jYo 



(28) 



From Eqs.||22J-(|2EJ> follows that in approximation 
the dynamics of inflationary cosmological models at a 
bounce does not depend on parameter a. In the case 
of models containing at a bounce ultrarelativistic matter 
and scalar fields, the dynamics in considered approxima- 
tion does not depend on ultrarelativistic matter, if inter- 
action terms in scalar field potential can be neglected. 

The analysis given in this Section shows that existence 
of limiting bound L and bounce surface in space P en- 
sures regular character of inflationary cosmological solu- 



tions in metrics, Hubble parameter and its time deriva- 
tive. Characteristic feature of inflationary cosmological 
models is the presence of the stage of regular transition 
from compression to expansion. The duration of this 
stage is several order smaller than duration of inflation- 
ary stage. If we take into account that duration of infla- 
tionary stage in chaotic inflation is extremely small Q], 
we can tell that our inflationary models correspond to 
regular Big Bang or Big Bounce. 



IV. PARTICULAR REGULAR INFLATIONARY 
COSMOLOGICAL MODEL 

As illustration of discussed theory regular cosmologi- 
cal models in the simplest particular case will be con- 
sidered in this Section. We will consider models includ- 
ing noninteracting scalar field with potential V(0) and 
ultrarelativistic matter (p m = ^p m )[3(|. In this case 

the bound L in space P is reduced to two cylindric sur- 

i 

faces = ± (4V + a -1 ) 2 • Bounce surface is reduced 
also to cylindric surfaces in the case under considera- 
tion, when the scale of extremely high energy densities 
is much smaller than the Planckian energy density (see 
(123). In connection with this we will consider instead 
of space P the plane of variables (cj), 0) and intersec- 
tions of the bound L and bounce surface with this plane. 
We have in this plane two bound L± -curves and in the 
case of flat models two bounce curves B\ and B>2 deter- 
mined by equation 39] 
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Each of two curves JSi^ contains two parts correspond- 
ing to vanishing of H+ or and denoting by (Bi + , 
B 2+ ) and B 2 ) respectively. If V is positive (neg- 

ative) in quadrants 1 and 4 (2 and 3) on the plane (</>, 
4>) 1 the bounce will take place in points of bounce curves 
B\+ and B 2 + (-Bi- and B 2 ) in quadrants 1 and 3 (2 
and 4) for i/ + -solutions (^--solutions) (see Fig. To 
analyze flat bouncing models we have to take into ac- 
count that besides regions lying between curves L± and 
corresponding bounce curves the sign of values H+ and 
H- for applying potentials is normal: H + > 0, H- < 0. 
The Hubble parameter H+ is negative in regions between 
curves (L + and (L- and B 2+ ), and the value of 

H_ is positive in regions between curves (L + and £>i_), 
(L_ and B 2 J). As it was noted above any cosmologi- 
cal solution has to contain both and i/ + -solution 
and regular transition from H- -solution to -H+-solution 
takes place in points of L± where H + = H-. In the 
case of open and closed models Eq. (J22J) determines 1- 
parametric family of bounce curves with parameter Rq. 
Bounce curves of closed models are situated in region be- 
tween two bounce curves B\ and B 2 of flat models, and 
in the case of open models bounce curves are situated in 
two regions between the curves: L + and B±, L_ and B 2 . 
In general case, when approximation ill'MI) is not valid, 
bounce surface in space P of cosmological models includ- 
ing scalar field and ultrarelativistic matter determined 
by Eq. (|20|l in space P depends on parameter a and it 
is not more cylindric surface. The situation concerning 
cosmological solutions of Eqs. ©-(EDI does not change. 

Below particular bouncing cosmological inflationary 
solution for flat model by using scalar field potential in 
the form V = \m 2 (f) 2 (m = 10~ 6 M p )is given. The solu- 
tion was obtained by numerical integration of Eqs. J7J), 
HI U| l and by choosing in accordance with Eq. H2Uf> the fol- 
lowing initial conditions at a bounce: 0o = 
(j) = V3.96757Vb {a = 10 14 A/" 4 ); radiation energy den- 
sity is negligibly small, initial value of i?o can be arbi- 
trary. The dynamics of the Hubble parameter and scalar 
field is presented for different stages of obtained bounc- 
ing solution in Figs. (by choosing M p — 1). The 
transition stage from compression to expansion (Fig. [2J 
is essentially asymmetric with respect to the point t = 
because of <j) ^ 0. In course of transition stage the Hub- 
ble parameter changes from maximum in module nega- 
tive value at the end of compression stage to maximum 
positive value at the beginning of expansion stage. The 
scalar field changes linearly, the derivative <j> grows at 
first being positive to maximum value <f> ~ <po and then 
the value of <f> decreases and becomes negative. In the 
case of presence of radiation, which can give the main 
contribution to energy density at a bounce, the duration 
of transition stage changes. Quasi-de-Sitter inflationary 
stage and quasi-de-Sitter compression stage are presented 
in Fig. [3J As was noted above, characteristics of infla- 
tionary stage do not depend practically on parameter a 



and coincide with that of GR. The amplitude and fre- 
quency of oscillating scalar field after inflation (Fig. 0J| 
are close to that of GR, however, the behaviour of the 
Hubble parameter after inflation in considering case with 
large value of parameter a is essentially noneinsteinian, 
at first the Hubble parameter oscillates near the value 
H = 0, and later the Hubble parameter becomes positive 
and decreases with the time like in GR. Before quasi-de- 
Sitter compression stage there are also oscillations of the 
Hubble parameter and scalar field not presented in Figs. 

Ultrarelativistic matter, which could dominate at a 
bounce has negligibly small energy densities at quasi-de 
Sitter stages. At the same time the gravitating matter 
could be at compression stage in more realistic bounc- 
ing models. As it follows from our consideration regular 
character of such inflationary cosmological models has to 
be ensured by cosmological equations of PGTG. 

The interaction between scalar fields and radiation 
leads to quantitative corrections of considered cosmolog- 
ical models. In accordance with Eq. (JHJ temperature 
corrections for scalar field potentials change the connec- 
tion between scale factor R and radiation energy density 
and can be also essential for more late stages of cosmolog- 
ical evolution, when energy densities are sufficiently small 
and consequences of PGTG and GR coincide. In partic- 
ular, these corrections can give additional contribution 
to the effect of acceleration of cosmological expansion (if 
relic scalar fields exist). 

Note that bouncing character have solutions not only 
in classical region, where scalar field potential, kinetic 
energy density of scalar field and energy density of grav- 
itating matter do not exceed the Planckian energy den- 
sity, but also in regions, where classical restrictions are 
not fulfilled and according to accepted opinion quantum 
gravitational effects can be essential. 

V. CONCLUSION 

As we see, the applying of gauge approach to gravita- 
tional interaction permits to build consequent field the- 
oretical scheme in the frame of 4-dimensional physical 
space-time, which is free of principal difficulties of GR by 
description of early Universe. Satisfying the correspon- 
dence principle with GR in the case of gravitating sys- 
tems with rather small energy densities, generalized cos- 
mological Friedmann equations lead to conclusion, that 
at extreme conditions gravitational interaction has the 
repulsion character in the case of usual gravitating sys- 
tems with positive values of energy density satisfying the 
energy dominance condition. This means, there is not 
necessity to refuse fundamental physical requirement of 
energy density positivity for physical matter. The solu- 
tion of PCS is obtained by classical description of gravi- 
tational field without quantum gravitational corrections. 
Moreover, from the point of view of developing approach, 



FIG. 1: Bound L±-curves and bounce curves for flat models in the case of potential V — ^m 2 (f> 2 . 




FIG. 2: The stage of transition from compression to expansion. 




FIG. 3: Quasi-de-Sitter stage of compression and inflationary stage. 



the Planckian era could be absent by evolution of our the transition stage from compression to expansion. To 
Universe. Unlike loop quantum cosmology, the Universe build a realistic cosmological model we have to know the 
is macro-object at all stages of its evolution including content and properties of gravitating matter at different 
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FIG. 4: The stage after inflation. 



stages of its evolution. It is of principal interest to inves- 
tigate physical processes at the beginning of cosmologi- 
cal expansion depending on limiting energy density and 
limiting temperature, and to obtain observable physical 
consequences depending on parameter a. From physical 
point of view, it is interesting also the building of cos- 
mological models by breakdown of its homogeneity and 
isotropy. 
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Because in the frame of gauge approach the gravita- 
tional interaction is connected with space-time trans- 
formations, the gauge treatment to gravitation has es- 
sential differences in comparison with Yang-Mills fields 
connected with internal symmetries groups. As a result, 
there are different gauge treatments of gravitational in- 
teraction not discussed in this paper. 
The structure of gravitational equations of PGTG for 
HICM does not depend on detailed form of quadratic 
part of gravitational Lagrangian Lq- As a result, these 
equations contain two independent parameters / and a, 
which are linear combinations of coefficients at terms of 
La quadratic in the curvature and torsion tensors. The 
mathematical requirement for cosmological equations for 
HICM to be differential equations of the same order as 
in GR can be fulfilled if a = that leads to Eqs (l)-(2). 
Unlike metric theories of gravitation, in PGTG terms 
of Lg quadratic in the curvature tensor do not lead to 
high derivatives in GCFE. Parameter a in Eqs. ©-©is 
defined as a = | (167r) 2 / MjT 4 . In our previous papers 
the parameter f3 — —a was used. 
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[35] The GCFE lead to limiting energy density also in the 
case a < if at extreme conditions p < |2^ |. 

[36] From mathematical point of view the appearance of spe- 
cific solutions is connected with the fact, that the Eqs. 
(l)-(2) were multiplied by Z 2 when its transformation 
to the (9)-(10) was performed. Specific solutions can be 
excluded by corresponding transformation of Eq. (10). 

[37] In GR a bounce is possible only in closed models if the 
following condition Vo — 4>o — h(PmO + 3p 7no ) > takes 



place. 

[38] If a -1 <C Mp the derivative H is negative in the neigh- 
bourhood of origin of coordinates in space P that leads 
to appearance of oscillating solutions of GCFE[3i1|. 

[39] The neighbourhood of origin of coordinates is not con- 
sidered in this approximation, the behavior of ex- 
tremum curves near origin of coordinates was examined 
in Ref. 



